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1 Introduction 

Cosmic ray (CR) particles arrive at the top of the Earth's atmosphere at 
a rate of around 10 3 per square meter per second. They are mostly ionized 
nuclei - about 90% protons, 9% alpha particles traces of heavier nuclei and 
approximately 1% electrons. 

CRs are characterized by their high energies: most cosmic rays are rela- 
tivistic, having kinetic energies comparable to or somewhat greater than their 
rest masses. A very few of them have ultra-relativistic energies extending be- 
yond 10 20 eV (tens of joules). 

In this series of lectures, delivered at the 2005 Mexican School of Astro- 
physics, an overview of the main experimental characteristics of the CR flux 
and their astrophysical significance is given. Particular emphasis is given to 
the upper end of the CR energy spectrum. 

Unfortunately, due to space limitations, only a fraction of the original 
content of the lectures is included in the present manuscript. In particular, 
the production mechanisms are not included and the fundamental topic of 
anisotropics is only dealt with in a very superficial way. 



2 Energy spectrum 

Thus, the cosmic ray energy spectrum extends, amazingly, for more than 
eleven orders of magnitude. All along this vast energy span, the spectrum 
follows a power law of index ~ 2.7. Therefore, the CR flux decreases approxi- 
mately 30 orders of magnitude from ~ 10 3 m 2 sec -1 at few GeV to ~ 1 km -2 
per century at 100 EeV. 

The only spectral features are a slight bending at around few PeV, known 
as the first knee, another at approximately 0.5 EeV known as the second knee, 
and a dip extending from roughly the second knee up to beyond 10 EeV, 



2 Gustavo Medina Tanco 

known as the ankle (see, figure^. Note that in the right panel the spectrum 
is multiplied by E 3 , a usual trick to highlight features that otherwise would 
be almost completely hidden by the rapidly falling flux. 



m -10 

I 

C 1 -1B 



Bdioo-i ana 
Satellite experiments 




£ 
> 



tfl 19 20 V\ 
log( ENERGY eV} 



KNEE 



2nd. ANKLE 




□ Akeno 1 km* 

* Tibet 

A Runjob 

* Proton Satell' 

* JACEE 



O AGASA 
A Haverah Park 
■Sr Yakutsk 
* Stereo Fly's Eye 

■ » - I 



1» it 20 21 

log(ENERGV in eV) 



Fig. 1. Cosmic ray energy spectrum and its main features: [left] a remarkably 
uniform power law with [right] few bends knee (few PeV), second knee ( 0.5 EeV), 
ankle (EeV to few tens of EeV) and the still poorly known highest energy tail. 
Adapted from pQ 



The second knee has been observed in the vicinity of 4 x 10 17 eV by Akeno 
|2, Fly's Eye stereo El El, Yakutsk H and HiRes 10 . The physical 
interpretation of this spectral feature is uncertain at present. It may be either 
the end of the Galactic cosmic ray component or the pile-up from pair creation 
processes due to proton interactions with the cosmic microwave background 
radiation during propagation in the intergalactic medium. 

The ankle, on the other hand, is a broader feature that has been observed 
by Fly's Eye -jSliniE] around 3 x 10 18 eV as well as by Haverah Park jS] at ap- 
proximately the same energy (3 x 10 18 eV). These results have been confirmed 
by Yakutsk [313] and HiRes [TUj. AGASA also observed the ankle, but they 
locate it at a higher energy, around 10 19 eV HJ. As with the ankle, more than 
one physical interpretations are possible, which are intimately related with the 
nature of the second knee. The ankle may be the transition point between the 
Galactic and extragalactic components, the result of pair creation by protons 
in the cosmic microwave background, or the result of diffusive propagation of 
extragalactic nuclei through cosmic magnetic fields. 
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Composition 

Certainly, one of the scientifically most relevant pieces of information inside 
the transition energy interval previously defined, is the precise chemical com- 
position of the primary CR flux as a function of energy. 

Several techniques have been used to determine the composition of cos- 
mic rays along the spectrum and, in particular, in the highest energy region 
|12| : (i) depth of maximum of the longitudinal distribution, X max |13II14| ; (ii) 
fluctuations of X max ^J^];(iii) muon density steepness of the 

lateral distribution function \IQ \T7\ ; (v) time profile of the signal, in particu- 
lar rise time of the signal |28j : (vi) curvature radius of the shower front; (vii) 
multi-parametric analysis, such as principal component analysis and neural 
networks |2H], etc.. Unfortunately, as is frequently the case in physics, when- 
ever several techniques are applied to measure the same physical magnitude, 
correspondingly, several results are obtained and, not always agreeable among 
themselves. As will be shown below, this is critical to the understanding of 
the astrophysics of ultra high energy cosmic rays. 

In order to analyze the astrophysical implications of the composition along 
the second knee/ankle energy region, it is more instructive to start from much 
lower energies. A significant point is the first knee. Following KASCADE |21j . 
a gradual change in composition is observed through the knee, from a lighter 
to a heavier composition. The first knee is a broad feature which can be 
understood as a composition of power law energy spectra with breaks that 
are in agreement with a rigidity scaling of the knee position. 

Therefore, at energies above few times 10 16 eV, the flux is dominated by 
iron nuclei. These particles are of Galactic origin and what is being detected 
is, very likely, the end of the efficiency of supernova remnant shock waves as 
accelerators as the Larmor radii or characteristic diffusion scale lengths of the 
nuclei become comparable to the curvature radius of the remnants, breaking 
down the diffusive approximation. If there are not more powerful accelerators 
in the Galaxy, the Galactic cosmic ray flux continue dominated exclusively 
by iron above 10 17 eV up to the highest energies produced inside the Milky 
Way. It must be noted that, even if the previous results are quantitatively 
dependent on the hadronic interaction model used in the data analysis, they 
are qualitatively solid and there is considerable degree of consensus on the ex- 
istence of a progressive transition in composition through the knee. At higher 
energies, the composition has been measured by several experiments in the 
past, e.g., Haverah Park, Yakutsk, Flys Eye, HiRes-MIA prototype and HiRes 
in stereo mode (see figure 

The X max data suggest that, above 10 166 eV, the composition changes 
once more progressively from heavy to light. At the lower limit of this energy 
interval, the composition is still heavy, i.e., iron dominated, in accordance to 
Kascade results. Nevertheless, at energies of 10 19 eV, even if still showing 
signs of a contamination by heavy elements, it is more consistent with a flux 
dominated by lighter elements. 
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Fig. 2. Variation of X max with energy (elongation rate) showing an apparent change 
in composition from heavy to light nuclei from ~ 30 PeV to ~ 30 EeV. This variation 
is at least possibly associated to the transition from Galactic to extragalactic cosmic 
rays. The lines indicate theoretical expectations corresponding to different hadronic 
interaction models. 



Despite the fact that there is a consensus among most of the experiments 
about the reality of this smooth transition, there is no consensus about the 
rate and extent to which the transition occurs. In fact, the combined data 
from the HiRes-MIA prototype and HiRes in stereo mode, signal to a much 
more rapid rapid transition from heavy to light composition (see figure yj , 
starting 10 17 eV but which would be over by 10 18 eV (301 ■ Beyond that point, 
the composition would remain light and constant. 

The later scenario, however, is not supported by the data of other experi- 
ments. Haverah Park, for example, shows a predominantly heavy composition 
up to 10 18 eV, followed by an abrupt transition to lighter values compati- 
ble with HiRes stereo at around 10 18 eV (see, figure 0Ja). Volcano Ranch, 
even though there is a single experimental point, is compatible with a heavy 
composition still at 10 18 eV, somewhat in accordance to Haverah Park data. 
Akeno (Al), on the other hand, is consistent with a continuation of the grad- 
ual transition from the second knee all across the ankle up to at least 10 19 
eV, only reaching there the same light composition that HiRes stereo claims 
from an order of magnitude below in energy. It must also be noted that above 
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10 19 eV AGASA is only able to set upper limits for the fraction of iron, but 
these limits are high enough to leave room for much more complex astrophys- 
ical scenarios with a substantial admixture of extragalactic ultra-high energy 
heavy elements |ol| . 

Figure ^b shows a compendium of several of the available measurements 
of composition between ~ 10 17 and ~ 10 19 eV, with their corresponding error 
bars, under the simplistic assumption of a binary mixture of protons and 
iron nuclei. The emerging picture is one complete uncertainty, which has deep 
practical implications and imposes severe limitations to theoretical efforts. 

At energies beyond ~ 10 EeV the composition is essentially unknown. 
However, it seems compatible with hadrons even if some photon contribution 
cannot be discarded [321 HUH EU (see figure (SJ. Although it is implicitly re- 
garded as purely protonic in many theoretical works, only upper limit exists 
for the iron fraction (e.g.,0Jb) and therefore not much can be said until more 
accurate measurements are made available by new generation experiments like 
Auger and TA. 
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Fig. 4. (a) Variation of the iron fraction inside the transition energy interval for var- 
ious experiments, (b) Idem, highlighting the uncertainties in composition inside the 
region encompassing the second knee and the ankle: literally almost any abundance 
of iron is allowed (adapted from \27\). 
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Fig. 5. Upper limits (95% CL) on cosmic-ray photon fraction for Auger |84j . 
AGASA (Al) [22, (A2) |S3 and Haverah Park (HP) 1351 EU data compared to 
some estimates based on TOP-DOWN models j^. (Reproduced from |84p. 



3 Galactic propagation 

The fundamental question of cosmic ray physics is, "Where do they come 
from?" and in particular, "How are they accelerated to such high energies?". 
These are difficult questions not fully answered after almost a century of 
history of the field. 
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Some very general hints can be obtained, however, through very simple ar- 
guments. The interstellar and intergalactic mediae are magnetized and, being 
charged, the CR are forced to interact with these fields. 

From the point of view of propagation of charged particles, the Galaxy 
behaves as a magnetized volume, where the field is structured on scales of 
kpc, with typical intensities of the order of some few micro Gauss. 

The Larmor radius of a nucleus of charge Ze can be conveniently parame- 
terized as: 



where EebV is the energy of the particle in units of 10 18 eV and rL,kpc is 
expressed in kpc. 

Equation Q clearly shows that, given the typical intensity of the mag- 
netic fields present in the interstellar medium (ISM), nuclei with energies 
below some few tens of EeV, regardless of their charge, have Larmor radii 
much smaller than the transversal dimensions of the magnetized Galactic 
disk. They must, therefore, propagate diffusively inside the ISM. This trans- 
forms the Galaxy in an efficient confinement region for charged particles with 
energies below the second knee. The confinement region is a flattened disk of 
approximately 20 kpc of radius and thickness of the order of a few hundreds 
of pc. 

Consequently, from the lowest energies and up to the second knee, the 
Galaxy is undoubtedly the source of the cosmic ray particles and of their 
kinetic energy. There is not a consensus about the actual source of the par- 
ticles in itself, but the two main lines of thought propose either nuclei pre- 
accelerated at the chromospheres of normal F and G type stars or ambient 
electrically charged nuclei condensed in the dense winds of blue or red giant 
stars [HB1 CHI • On the other hand, several acceleration mechanism must be at 
play but it is widely expected that the dominant one is first order Fermi accel- 
eration at the vicinity of supernova remnant shock waves. Nevertheless, the- 
oretically, the Galactic accelerators should become inefficient between ~ 10 17 
to ~ 10 18 eV. This upper limit could be extended to ~ 10 19 eV if additional 
mechanisms were operating in the Galaxy, e.g., spinning inductors associated 
with compact objects or cataclysmic events like acceleration of iron nuclei 
by young strongly magnetized neutron stars through relativistic MHD winds 



At energies above the second knee, particles start to be able to travel from 
the nearest extragalactic sources in less than a Hubble time. Consequently, 
at some point above 10 17 ' 5 eV a sizable cosmic ray extragalactic component 
should be detectable and become dominant above 10 19 eV. Therefore, it is 
expected that the cosmic ray flux detected between the second knee and the 
ankle of the spectrum be a mixture of a Galactic and an extragalactic flux, 
highlighting the astrophysical richness and complexity of the region. 




(1) 
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The type of propagation strongly depends on the charge of the correspond- 
ing nucleus. 

Protons with energies > 10 17 eV have gyroradii comparable or larger than 
the transversal dimensions of the effective confinement region and, therefore, 
can easily escape from the Galaxy. On the other end of the mass spectrum, 
just the opposite occurs for iron nuclei that, even at energies of the order of 
10 19 eV, have gyroradii < 10 2 pc and must be effectively confined inside the 
magnetized interstellar medium. 

The previous results are based only on the consideration of the regular 
component of the Galactic magnetic field. However, there exist a superimposed 
turbulent component whose intensity is at least comparable to that of the 
regular field. Its spectrum seems to of the Kolmogorov type, extending from 
the smallest scaled probed, ~ 10° pc, to L c ~ 100 pc, the correlation length 
of the turbulent field. 

Wave-particle interactions between cosmic rays and MHD turbulence are 
resonant for wavelengths of the order of the Larmor radius, A ~ ri,. This 
means that, for a nucleus of charge Z, a critical energy can be defined, 



bellow which modes resonant with the particle gyroradius exist that are able 
to efficiently scatter the particle in pitch angle. Consequently, at energies 
bellow E c the diffusion coefficient is small enough for the particles trajectory 
to be diffusive. At energies above E c , on the other hand, the propagation is 
essentially ballistic. 

Due to the interaction with the turbulent magnetic component, protons 
experience a propagation regime very different than iron nuclei inside the ankle 
energy region. Protons propagate ballistically in the interstellar medium above 
~3x 10 17 eV, while iron nuclei propagate diffusively even at energies > 10 19 



Therefore, along the energy region extending from the second knee up to 
almost the end of the ankle, all nuclei from p to Fe, i.e. 1 < Z < 26, experi- 
ence a transition in their propagation regime inside the interstellar medium 
changing gradually from diffusive to ballistic as the energy increases. 

A pictorial example of how this transition takes place can be seen in figure 
IHla-d |41| . where it is shown how protons with energies ranging from 0.5 to 6 
EeV injected at the Galactic center propagate out of the Galaxy assuming a 
characteristic BSS magnetic model. 

It must be noted that, despite the fact that the deflections induced by the 
Galactic magnetic field (GMF) diminish rapidly with energy for all nuclei, they 
can still be important even at the highest energies. This is more critical for 
heavier primaries and for all nuclei traversing the central regions of the Galaxy. 
Figurc[7|illustratcd this point by showing the intrinsic deflections experienced 
by proton and iron nuclei as a function of arrival Galactic coordinates |42| . 




L c w 10V =>■ E c ,EeV ~ 0.5 x Z (2) 
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Fig. 6. Changes in propagation regime inside the Galaxy at energies of the second 
knee and ankle. 



Figure [3] shows Galactic deflections for 100 EeV protons for a more sophis- 
ticated GMF model inspired on Hans proposal . The GMF is modelled by 
a disk BSSS component of 100 pc half thickness, embedded in an ASSA halo 
and a dipolar component originated in a Southward magnetic momentum an- 
chored at the Galactic center. Galactic coordinates in an Aitof projection are 
used. 
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Fig. 7. Intrinsic deflections due to the GMF (BSSS model) suffered by protons 
and Fe nuclei at 4 x 10 19 eV and protons at 2.5 x 10 20 eV as a function of arrival 
direction. Galactic coordinates are used. 



4 Extragalactic propagation 

4.1 Superposition of the extragalactic and Galactic fluxes 

In the same way as the magnetic characteristic of the interstellar medium allow 
Galactic particles at these energies to escape into the extragalactic environ- 
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ment, extragalactic cosmic rays are also able to penetrate inside the Galactic 
confinement region. But, of course, extragalactic particles must first be able 
to reach us from the nearest Galaxies in less than a Hubble time. 

A crude approximation to this effect can be made in the following way. 

Faraday rotation measurements statistically impose to the extragalactic 
magnetic field the following restriction |43| : 

B x h]J 2 < InG x Mpc 1/2 (3) 

where L c is the correlation length of the magnetic field that we assume, some- 
what arbitrarily, as being of the order of 1 Mpc. Assuming that the diffusion 
coefficient can be estimated by the Bohm approximation, 

K « ir L c (4) 

and using equation Q the diffusion coefficient can be written: 

0.1 (E EeV \ Mpc 2 
KK ^y-B^) My7 (5) 

The diffusive propagation time from an extragalactic source at a distance 
D can be estimated as: 



D 2 

T W 

K 

or, using equation (J5J: 



(6) 



r„„, « H) :D 2 Mpc xZx ( J^) (7) 
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Equation Q shows that there is a rather restrictive magnetic horizon. Basi- 
cally, no nucleus with energy smaller than 10 17 eV is able to arrive from regions 
external to the local group [D ~ 3 Mpc). Taking as a minimum character- 
istic distance D — 10 Mpc, which defines a very localized region completely 
internal to the supergalactic plane and even smaller than the distance to the 
nearby Virgo cluster, only protons with E > 2 x 10 17 eV, or Fe nuclei with 
E > 5 x 10 18 eV are able to reach the Galaxy in less than a Hubble time. 

Therefore, it is at the energies of the second knee and the ankle that 
different nuclei start to arrive from the local universe. Concomitantly, at these 
same energies, the magnetic shielding of the Galaxy becomes permeable to 
these nuclei, allowing them to get into the interstellar medium and, eventually, 
to reach the solar system. Effectively, the energy interval from ~ 2 x 10 17 
to 10 19 eV is the region of mixing between the Galactic and extragalactic 
components of cosmic rays. 

Above few times 10 17 eV, the dominant interactions experienced by cosmic 
rays are due to the cosmic microwave background radiation (CMBR) and, 
additionally in the case of nuclei, to the infrared background (CIBR) |44j . 
The diffuse background in radio, despite its much lower density, must in turn 
become important at high enough energies. 

At energies above ~ 10 192 eV, the dominant process is the photo produc- 
tion of pions in interactions with the CMBR (see figure^a), which drastically 
reduces the mean free path of protons to e few Mpc, making the universe op- 
tically thick to to ultra-high energy cosmic rays ijHlb). This interaction, in the 
most conservative models, should produce a strong depression in the energy 
spectrum, with a major fall in the observed flux above 10 20 eV, the so called 

gzk cut-off mum. 

At energies smaller than ~ 10 19 ' 2 eV, the dominant process is the photo- 
production of electron-positron pairs in interactions with the CMBR. At these 
lower energies the attenuation length attain values of the order of Gpc and, 
therefore, the universe is essentially optically thin to energetic baryons. CR 
observations at these energies, sample the universe at cosmological distances, 
contrary to the highest energies, that only sample a sphere of a few tens of 
Mpc in diameter, a small portion of the local universe |5()l 1511 152] . Therefore, 
strictly from the point of view of propagation in the extragalactic medium, in 
going down from the highest energies to the transition region, the observable 
horizon drastically increases from 10 1 to 10 3 Mpc, i.e., essentially the whole 
universe. 

It can also be seen from figure a that the dependence of the cross section 
with energy is suggestive, since its shape resembles that of the ankle in the 
cosmic ray energy spectrum. In fact, the structure of the ankle can be ex- 
plained exclusively as a result of pair photo-production by nucleons travelling 
cosmological distances between the source and the observer |4*5], 

The energy region where the superposition of the Galactic and extragalac- 
tic spectra takes place is a theoretically challenging region, where the smooth 
matching of the two rapidly varying spectra has yet to be explained. It must 
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Fig. 9. (a) Cross section for pair production and pion production in interaction 
with the CMBR. The positions of the second knee and the ankle are also shown, 
demonstrating that electron positron pair production is the relevant interaction in 
the Galactic-extragalactic transition region. Note the similitude between the shape 
of the cross section for this interaction and the shape and location of the ankle, (b) 
Attenuation length in Mpc as a function of energy |45|. showing how the universe, 
which is opaque for at energies above the photo-pion production threshold, becomes 
transparent to lower energy barions. 



be noted that, even if the shape of the spectrum is important, it is by far 
insufficient to decipher the the underlying astrophysical model. The Galactic 
magnetic fields are intense enough to dilute any directional information, which 
prevents the discrimination among the galactic and extragalactic components 
from the arrival direction of the incoming particles. The variation of the com- 
position as a function of energy turns then into the key to discriminate both 
fluxes and to select among a variety of theoretical options. 

As in the case of the interstellar medium, it is expected that the intergalac- 
tic medium has a strong magnetic turbulent component which can severely 
affect propagation 02J EHI HSU EH EH| • The correlation length estimated 
from Faraday rotation measurements, L C1 is consistent with a maximum wave- 
length for the MHD turbulence determined by the largest kinetic energy in- 
jection scales in the intergalactic medium, L max ~ L c ~ 1 Mpc. Therefore, 
analogous to equation J5J: 



z 



EeV 



nG 



IMpc E CtEeV wl.OxZ 

(8) 



14 



Gustavo Medina Tanco 




♦ Agasa 



Energy [eV] 



Fig. 10. Total cosmic ray spectrum Irom the combined data of several experiments. 
From the theoretical point of view, the transition region is highly complex and the 
Galactic and extragalactic models suffer their most critical test since the correspond- 
ing fluxes must be simultaneously matched both in intensity and energy. (Adapted 
from: 1351 .1 



For a given nucleus of charge Ze, the propagation is ballistic for E > E c 
being diffusive otherwise. Therefore, protons are ballistic above ~ 10 18 eV, 
but diffusive at the energies of the second knee. Iron nuclei, on the other 
hand, propagate diffusively along the ankle and even at energies as high as 
~ 5 x 10 19 eV. The boundaries for the transition between the ballistic and 
diffusive propagation regime for proton and Fe nuclei are shown in the figure 

EH 

Furthermore, besides the total intensity and the minimum wavenumber, 
also the energy distribution among the different modes, that is the type of 
turbulence present in the intergalactic medium, have observational expression. 
In this case, the affected portion of the extragalactic spectrum is the region 
of lower energies, where the flux is strongly suppressed by magnetic horizon 
effects. Figure 1771 shows clearly this effect for three different assumptions for 
the diffusion coefficient. 

Obviously, this has profound theoretical implications not only for the struc- 
ture of the extragalactic magnetized medium, but also for cosmic ray accel- 
eration conditions inside the Galaxy. This is exemplified in figure ^] where it 
is graphically illustrated that, by subtracting a given extragalactic spectrum 
from the observed total spectrum, conclusions can be drawn about relevant 
aspects of the Galactic component. For example, an extragalactic spectrum 
that has a small contribution at low energies, can imply the existence of ad- 
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Fig. 11. Correlation between the detailed structure of the lower end of the extra- 
galactic spectrum and the type of turbulence present in the intergalactic medium. 
This emphasize the importance of the intergalactic turbulent component in the ob- 
served matching between the Galactic and extragalactic cosmic ray flux. (Adapted 
from: |55|.1 



ditional acceleration mechanisms in the Galaxy other than the shock waves 
of supernova remnants. 

4.2 The highest energies 
Propagation of protons 

As was mentioned in the previous section ETTl fsee figureEJ, above the threshold 
for photo-pion production by protons interactiing with the CMBR, ~ 40 EeV 
the universe becomes rapidly opaque for hadrons as the attenuation length 
goes down to values as low as ~ 10 Mpc at few x 10 20 eV. This determines a 
relatively small (in cosmological terms) maximum distance scale, Rgzk — 50- 
100 Mpc, to the sources that are able to contribute appreciably to the detected 
CR flux. 

Under very general assumptions regarding the nature of the primaries and 
the cosmological distribution of the sources, photo-pion production should 
lead to the formation of a pile-up immediately followed by a severe reduction 
in CR flux, popularly known as the GZK cut-off. The existence of this spectral 
feature was proposed short time after the discovery of the CMBR 021 EH1 but 
its actual existence is still a matter of considerable debate. 
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Fig. 12. Impact of the detailed characteristics of the extragalactic spectrum on 
our comprehension of the most powerful acceleration mechanisms in our Galaxy. 
(Adapted from: gO].) 



At present, there are conflicting measurements coming from two different 
experiments: AGASA and HiRcs. The first one is a surface detector while the 
second one is a fluorescence detector, which further complicates the compari- 
son of their results. 

As shown in figure lT3l the differences are not only quantitative but, funda- 
mentally, qualitative. While HiRes apparently shows the expected GZK flux 
suppression, AGASA seems incompatible with this result, showing an energy 
spectrum that extends undisturbed well beyond 100 EeV. There is also an ap- 
parently large difference in flux between both experiments at energies below 
the cut-off. However, the fact that the energy spectrum has been multiplied by 
E 3 in figure 1131 should be taken into account when assessing the significance 
of such difference. 

Actually, both results might be reconciled at the 1.5 a level by re-scaling 
the energy of either of the experiments by 30% or by 15 % each The 
Auger Observatory, being the largest detector ever built and having hybrid 
capacity (i.e., simultaneous fluorescence and surface detection |fc>3p has the 
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Fig. 13. Comparison between the AGASA and Hires monocular spectra. (Adapted 
from: |rJT].) 



potential, but not enough statistics at present, to give a definite answer to 
this fundamental problem |64| . 

The absence of the GZK cut-off, if confirmed, could be compatible with 
a wide range of astrophysical scenarios. At least three possibilities can be 
considered, some rather conservative, some more exotic: 

• The distance scale between sources could be large enough that, by chance, 
the few (or single) sources inside the GZK sphere dominate the flux, the 
rest of the population of UHECR accelerators being too distant to con- 
tribute appreciably to the observed flux at Earth. 

• The primary CR might be particles that do not interact with photons 
or do so at much larger, and as yet unobserved, energies. These could 
be familiar standard model particles that present unexpected behavior at 
ultra-high energies, like neutrinos with hadronic cross section that can de- 
velop showers in the atmosphere resembling those expected from proton 
primaries 65,66. Another possibility could be a new stable hadron, heav- 
ier than a nucleon, for which the threshold for photo-meson production 
would be at higher energy. An example of the latter would be uhecrons, 
e.g., a uds-gluino bound state |67| . 

• The primaries might be normal hadrons, but Lorentz invariance, which 
has never been previously tested at 7 ~ 10 11 , could be violated at ultra- 
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high energies, hampering photo-meson production [681 1691 170] . The small 
violation of Lorentz invariance required might be the result of Planck scale 
effects IHIIZll. 

• The observed spectrum could be the superposition of two components: (a) 
a hadronic component with a GZK cut-off and (b) a harder, top-down 
component that becomes dominant above ~ 100 EeV. These second com- 
ponent could be originated in the decay or annihilation of super heavy 
dark matter or topological defects [23 EJ ESI E] ■ These scenarios have 
the general disadvantage of overproducing ultra-high energy neutrinos and 
photons rather than nucleons, which seems to be increasingly constrained 
by the observation Nevertheless, there could still be models, like 

those involving necklaces, that could present an appreciable baryon con- 
tent at energies > 100 EeV 76 . In any case, these models suffer from a 
discomforting level of fine tuning with respect to the normalization of the 
intensities of the GZKed and the top-down spectra. 

It must also be noted that the presence of the GZK flux suppression does 
not imply the non-existence of supra-GZK particles in the CR flux. These 
have certainly been detected by at least Volcano Ranch |781 179) . Flys Eye 
|80j . AGASA, HiRes and, more recently, Auger [SJ. This means that, the 
detection of the GZK feature does not solve the puzzle about the generation 
of UHECR. 



Propagation of UHE photons 

The propagation of photons is dominated by their interaction with the photon 
background. The main processes are photon absorption by pair-production on 
background photons (77b — > e + e~), and inverse Compton scattering of the 
resultant electrons on the background photons. These two processes acting in a 
chain are responsible for the rapid development of electromagnetic cascades in 
the intergalactic or interstellar media, draining energy to the sub-TeV region. 

For a given UHE photon of energy E 11 the minimum background photon 
energy, E^, for electron-positron pair production is: 

ml 2.6 x 10 11 
E b = ~ — eV (9) 

and the corresponding cross section peaks near the threshold: opp cx (m^/s)* 
ln(s/2m e ). Inverse Compton scattering, on the other hand, has no threshold 
but its cross section is also largest near the 77b pair production threshold. 
Therefore, the most efficient background for both processes is given by equa- 
tion Jj|J). For UHE this means that the cosmic radio background, whose mag- 
nitude is highly uncertain, is dominant followed by the CMBR below 10 17 - 
10 18 eV. At progressively lower energies, the CIRB and optical background 
are important. 
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In the Klein-Nishina limit, s m^, one of the components of the 77b - 
produced pair carries most of the energy of the energy of the UHE photon. 
This leading particle, afterwards, undergoes Compton scattering in the same 
limit, for which the inelasticity is very near 1. Therefore, the Compton up- 
scattered photon still has an appreciable fraction of the energy of the original 
UHE photon. The presence of magnetic fields in the medium may speed up the 
development of the cascade by draining the electron and positron energy due 
to synchrotron radiation. The larger the fields, the smaller the penetration. 

The electromagnetic cascades produced in this way can propagate an ef- 
fective distance that is much larger than the interaction length yet, severely 
limit our UHE-7 horizon to the nearest regions of the supergalactic plane. Fig- 
ure El 73] shows the effective penetration length of electromagnetic cascades 
for two different estimates of the cosmic radio background and two different 
average intergalactic magnetic field intensities. 

Since the single pair cross section decreases as ln(s)/s for s 3> to 2 , multi- 
ple pair production becomes important at extreme energies. Thus, double pair 
production (77^ — > e+e~e+e~) begins to dominate above ~ 10 21 - 10 23 eV. 
The relevant process for electrons is triple pair production (e ± 7(, — > e ± e _ e + ), 
whose attenuation becomes dominant at ~ 10 22 eV. Other processes (e.g., 
moun, tau or pion pair production, double Compton scattering, gamma scat- 
tering and pair production of single photons in magnetic fields) are in general 
negligible for electromagnetic cascade development. However, at energies in 
excess of 10 24 eV, the pair production of single photons in the Galactic mag- 
netic field should eliminate all the photons above that energy from specific 
lines of sight, generating an arrival direction anisotropy. 

The penetration lengths shown in figure PHI combined with the threshold 
energies given by equation (JjJJ, imply that the injection of UHE-7 in the in- 
tergalactic medium results in the pile-up of photons at energies below 100 
MeV, whose contribution to the diffuse cosmic 7 background is already obser- 
vationally constrained by EGRET. This overproduction of low energy diffuse 
photons is a strong restriction for top-down UHECR production models. 

Besides the limitations imposed by the possible overproduction of low en- 
ergy photons, the results in figure HH have other profound implications for top- 
down scenarios. Models that claim decay or annihilation of DM, in general, 
tend to produce mainly photons and only e few percent baryons. Therefore, 
in such models, most of the detected CR should be photons from our own 
Galactic halo, with perhaps some localized contribution from Andromeda. 
The products of the decay of more distant dark matter would be cleared 
from photons, and only the small fraction of remaining baryons, suppressed 
by GZK effects, would give a positive contribution at Earth. In any case, for 
most top-down models, photons should be dominant above the GZK cut-off, 
but with perhaps a sizable baryonic component. 
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Fig. 14. Penetration of electromagnetic cascades in the intergalactic medium. Solid 
lines correspond to a fiduciary value of 1 nG for the IGMF and dotted lines to a 
very low IGMF, 10~ nG. Thick and thin lines correspond to different estimates of 
the cosmic radio background. (Adapted from: 



Propagation of nuclei 

Heavy nuclei are attenuated by basically to processes: photodisintegration and 
electron-positron pair creation b y interaction with background photons. 

For the same total energy, the threshold photo-pion production for a nu- 
cleus of mass A increases to E t h — 4 x 10 19 x A. Therefore, given the energies 
observed at present, pion production is not relevant for nuclei heavier than 
He. 

Figure El shows that below ~ 20 EeV, all nuclei are able to travel for 
virtually a Hubble time, while Fe can do the same up to ~ 100 EeV. Above 
that energy, nuclei start to disintegrate fast and the loss time are highly 
reduced. 

It is clear from composition observations around the first knee of the cosmic 
rays spectrum that the acceleration mechanism by shock waves, either first 
order Fermi or drift are limited by the magnitude of the radius of curvature 
of the shock. This results in the preferential acceleration of large charge (eZ) 
nuclei, those with the smallest Larmor radius, to the highest energies. 

Even if the mechanism responsible for the acceleration of the ultra-high 
energy extragalactic component is essentially unknown, the most conservative 
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Fig. 15. Energy loss time (right axis: length) vs. energy for photodisintegration 
on background photons: radio, CMB and IR. Helium, Carbon, Silicon and Iron are 
shown. Single, double and multi-nucleon emissions are included. (Adapted from: 

E3-) 



view points to bottom-up mechanisms. If the later is actually the case, then 
the most economic assumption is that chock wave acceleration. 

In this minimalistic, but still realistic, scenario the most likely high energy 
output would be heavy nuclei, very likely Fe as in the Galactic case. At lower 
energies, progressively lighter nuclei should be observed due to two factors: (a) 
the acceleration process in itself and (b) the photo-disintegration on flight of 
the heavier nuclei due to their interaction with the CIBR. The latter process 
is very efficient, and can extract approximately one nucleon every few Mpc at 
the highest energies, depending on the CIBR level. Since photodisintegration 
occurs, to a good approximation, at constant energy per nucleon, disintegra- 
tion of a nucleus A at energy E will produce light nuclei at energies nE/A 
and (A — n)E/A with preferentially small n (e.g., n = 1, 2, 3, 4). 

Figure ITrjl shows that power law spectra injected at cosmological sources 
with different compositions can produce experimentally very similar at the 
highest energies. Nevertheless, they can always be distinguished at smaller 
energies in the ankle region. 

In particular, figure [TBI a showshow a purely protonic flux can reproduce 
the ankle feature solely as an effect of photo production of electron-positron 
pairs in interactions with the CMBR. In this case, the transition between the 
Galactic and extragalactic fluxes must be located at the second knee or very 
near to it. 
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Fig. 16. Different primary compositions may produce extragalactic spectra essen- 
tially indistinguishable at the highest energies. However, as shown in (a) for protons 
and (b) for heavier mix, the spectra are considerably different at lower energies 
below the bottom of the ankle. (Adapted from: [84].) 



Figure figure ^j]b, on the other hand, shows that, for a heavier mixed 
composition, the extragalactic spectrum falls down steadily with decreasing 
energy. In this scenario, the ankle must be the result of the composition be- 
tween the Galactic and extragalactic spectra. Moreover, the composition will 
be a strong function of energy inside this interval, giving an additional tool 
to assess details of the astrophysical model. 

The effects of photo-disintegration are clearly shown in figure |85| , where 
the evolution of a pure iron injected spectrum as it propagates out from the 
source. Blue histograms correspond to p, white to the original surviving Fe 
and red to intermediate mass nuclei. 

As the distance to the source increases, intermediate nuclei are produced 
at increasingly smaller masses and less total energy (fragmentation takes place 
at roughly constant mass per nucleon) . At distances of the order of the GZK 
horizon, the region with larger mixing of nuclei , i.e., with a larger composition 
gradient, is the ankle. Consequently, this is the ideal region for the discrimi- 
nation of the primary composition from local composition measurements. 



5 Cosmic magnetic fields and anisotropy 

Luminous matter, as traced by galaxies, as well as dark matter, as traced by 
galaxies and clusters large scale velocity fields, is distributed inhomogeneously 
in the universe. Groups, clusters, superclusters, walls, filaments and voids are 
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Fig. 17. Variation of the composition as a function of distance to the source for 
pure Fe power law injection. Blue histograms correspond to p, white to the original 
surviving Fe and red to intermediate mass nuclei. (Adapted from: |85| .) 



known to exist at all observed distances and are very well mapped in the local 
universe. Hence, the distribution of matter inside the GZK-sphere is highly 
inhomogeneous and so is, very likely, the distribution of UHECR sources. 

Synchrotron emission and multi-wavelength radio polarization measure- 
ments show that magnetic fields arc widespread in the Universe. But how do 
they encompass the structure seen in the distribution of matter we do not yet 
know 

The available limits on the IGMF come from rotation measurements in 
clusters of galaxies and suggest that Bxgm x L^J 2 < 1CT 9 G x Mpc 1 / 2 
where L c is the field reversal scale. Note, however, that this kind of mea- 
surement doesn't set an actual limit to the intensity of the magnetic field 
unless the reversal scale is known along a particular line of sight. The latter 
means that, depending on the structure of the IGMF, substantially different 
scenarios can be envisioned that are able satisfy the rotation measurement 
constraints. 

Unfortunately, we do not know what is the actual large scale structure 
of the IGMF. Nevertheless, we can imagine two extreme scenarios that are 
likely to bound the true IGMF structure. In figure ^] calculations of large 
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Fig. 18. two possible extreme models for the IGMF structure: (top) laminar struc- 
ture and (bottom) cellular structure. These schematic plots are adaptations by hand 
made on top of IGMF and density calculations by |1UU| 

scale structure formation by Ryu and co-workers 100 have been modified by 
hand to exemplify these scenarios. The top frame displays Ryu's IGMF sim- 
ulation results in the background showing how, by z = 0, the magnetic field 
has been convected together with the accretion flows into walls, filaments and 
clusters, depleting the voids from field. According to these calculations, the 
magnetic field is confined in high density, small filling factor regions, bounded 
by a rather thin skin of rapidly decreasing intensity, surrounded by large vol- 
umes of negligible IGMF. As suggested by the free-hand lines on top of the 
figure, the IGMF inside structures is highly correlated in scales of up to tens 
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of Mpc. Furthermore, in order to comply with the rotation measurement con- 
straints mentioned before, the intensity of the magnetic field inside the density 
structures must be correspondingly high, 0.1-1 fiG, which is comparable with 
GMF values within the interstellar medium. We will call the latter scenario 
laminar- structure . 

The second model, that we will call cellular- structure , is depicted in the 
bottom panel of figure El We imagine the space divide into adjacent cells, 
each one with a uniform magnetic field randomly oriented. We identify the 
size of a cell with the magnitude of the local reversal scale. Furthermore, one 
can assume that the intensity of the magnetic field scales as some power of the 
local matter (electron) density and, consequently, the rotation measurement 
constraint Bigm x lI^ 2 < 10~ 9 G x Mpc 1 / 2 tells how the reversal scale, 
i.e., the size of the cells, should be scaled. A convenient reference, such as the 
IGMF in the Virgo 86 or Comma jHOJ cluster can be used for normalization. 
The cellular-structure scenario leads to a more widespread IGMF, filling even 
the voids. The observational constraints imply then that the IGMF varies 
much more smoothly, from 10 _10 G inside voids to a few times 10 -9 -10 -8 G 
inside walls and filaments, only reaching high values, 0.1 — 1 /iG, inside and 
around clusters of galaxies. 

Observations are not enough at present to distinguish between these two 
scenarios. Nevertheless, we can still try to asses what are their implications 
for UHECR propagation which, inspecting figure IT§1 must be important. 



Fig. 19. Typical Larmor radii of nuclei in both IGMF scenarios, showing the very 
different scales involved in each model. 
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Fig. 20. (a) Simplified model of a wall, or slab, containing UHECR sources and 
surrounded by a void. The magnetic field configuration is representative of the lam- 
inar model, (b) Cross section of the wall in (a) at the plane z = 0. Several particles 
trajectories are shown for proton injection at E = 100 EeV. Adapted from 94 . 

5.1 UHECR propagation in a laminar-IGMF 

This is the most difficult scenario to dealt with because it does not accept 
a statistical treatment and results are very dependent on details about the 
exact magnetic field configuration inside the GZK-sphere, which is beyond 
our present knowledge. 

A simpler approach is to study the UHECR emissivity of a single wall 
surrounded by a void I94j. Figure 1201 a shows the corresponding model 
for a wall immersed in a void; the magnetic field inside the wall has two 
components, o uniform field along the z-axis of intensity 0.1 /iG and a random 
component with a Kolmogoroff power law spectrum of amplitude equal to 30% 
of the regular component. One hundred UHECR sources are included inside 
the wall, and each one of them injects protons at the same rate and with 
the same power low energy spectrum, dNi n j/dE oc E~ 2 . Pair-production and 
photo-pion production losses in interactions with the infrared and microwave 
backgrounds respectively are also included. The wall has a radius of 20 Mpc, a 
thickness of 5 Mpc and is sandwiched by a transition layer 5 Mpc in thickness 
where the magnetic field decreases exponentially up to negligible values inside 
the surrounding void. Once the system reaches steady state, a detector can 
be shifted around the wall to simulate observers at arbitrary positions with 
respect to the wall. In a real situation, this system could be representative, 
for example, of the supergalactic plane; in that case the Milky Way, i.e. we, 
the observers, should be located at some point on the x-z plane (but we don't 
know at what angle with respect to the z-axis). The simulations show that 
the UHECR flux measured can vary by three orders of magnitude depending 
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of the relative orientation between the wall, the field and the observer. At the 
same time, almost all directional information is lost, and the strength of the 
GZK-cut-off would vary considerably as a function of orientation . 

The previous effects can be intuitively understood by looking at figure 
1201b. which shows a cross section of the wall in figure |2~01 a at the plane z = 0. 
Several particles trajectories are shown for proton injection at E = 100 EeV, 
with different azimuthal angles and a slight elevation with respect to the x- 
y plane. It can be seen that there is nothing like a random walk: particles 
tend to be trapped inside the wall and move in a systematic way. Most of the 
particles drift perpendicularly to the regular field while their guiding centers 
bounce along the field. It can also be seen how the gyroradii decrease as the 
particles lose energy in interactions with the radiation backgrounds. Even the 
few particles that do escape from the wall, do so in an non-isotropic manner 
(e.g., predominantly to the right for y > 0). 

The laminar IGMF model is, actually, the worst scenario for doing some 
kind of astronomy with UHECR. It would be very difficult to interpret the 
UHECR angular data and to identify individual particle sources. Furthermore, 
the significance of any statistical analysis would be greatly impaired due to 
systematics. 

Further studies on this model can also be found, for example, in references 

una una. 

5.2 UHECR propagation in a cellular-IGMF 

The cellular model is the easiest scenario to deal with numerically and, by 
far, the most promising from the point of view of the astrophysics of UHECR. 
This is also the IGMF model that has been used probably more frequently in 
the literature EZl CU! IMl 0121 0E3 EH EH EH] ■ 

The main assumption here, is that the intensity of the magnetic field scales 
with density. Indeed, for those spatial scales where measurements are avail- 
able, the intensity of the magnetic fields seem to correlate remarkably well 
with the density of thermal gas in the medium. This is valid at least at galac- 
tic and smaller scales |101l Wf\ . It is apparent that B can be reasonably well 
fitted by a single power law over ~ 14 orders of magnitude in thermal gas 
density at sub-galactic scales. A power law correlation, though with a differ- 
ent power law index, is also suggested at very large scales (c.f. figure 5 in 
Medina Tanco 2000b), from galactic halos to the environments outside galaxy 
clusters, over ~ 4 orders of magnitude in thermal gas density. This view |101j 
is, however, still controversial 02] ■ I n ^ct, magnetic fields in galaxy clusters 
are roughly ~ 1 /iG, which is of the order of interstellar magnetic fields; fur- 
thermore, supracluster emission around the Coma cluster suggests /iG fields 
in extended regions beyond cluster cores. The latter could indicate that the 
IGMF cares little about the density of the associated thermal gas density, 
having everywhere an intensity close to the microwave background-equivalent 
magnetic field strength, Bbge — 3 x 10~ 6 G. 
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Taking the view that a power law scaling exist, a model can be devised 
in which the IGMF correlates with the distribution of matter as traced, for 
example, by the distribution of galaxies. A high degree of non-homogeneity 
should then be expected, with relatively high values of Biqmf over small 
regions (< 1 Mpc) of high matter density. These systems should be immersed 
in vast low density/low Bigmf regions with Bigmf < 1CP 9 G. Furthermore, 
in accordance with rotation measurements, the topology of the field should 
be such that it is structured coherently on scales of the order of the correla- 
tion length L c which, in turn, scales with IGMF intensity: L c oc Bj^ MF (r). 
Bigmf should be independently oriented at distances > L c . Therefore, a 3D 
ensemble of cells can be constructed, with cell size given by the correlation 
length, L c , and such that: L c oc Bj GMF (r), while Bigmf oc Pgai 5 { r ) (K)lj or 

oc p 2 Jai( r ) (f° r frozen-in field compression), where p ga i is the galaxy density, 
and the IGMF is uniform inside cells of size L c and randomly oriented with 
respect to adjacent cells The observed IGMF value at some given 

point, like the Virgo cluster, can be used as the normalization condition for 
the magnetic field intensity. The density of galaxies, p ga i, is estimated using 
cither rcdshift catalogs [like the CfA Redshift Catalogue (nHUMj or the PSCz 
|8ti|]. or large scale structure formation simulations The latter is a con- 
venient way to cope with, or at least to assess the importance of, the several 
biases involved in the use of galaxy redshift surveys to sample the true spatial 
distribution of matter in 3D space. 

The relevant energy losses for UHECR during propagation are: pair pro- 
duction via 7-7 with CMB for photons, redshift, pair production and photo- 
pion production in interactions with the CMB for nuclei and, for heavy nuclei, 
also photo-disintegration in interactions with the infrared background. All of 
these can be appropriately included [SSI 11031 [73] . 

The spatial distribution of the sources of UHECR is tightly linked to 
the nature of the main particle acceleration/production mechanism involved. 
However, in most models, particles will either be accelerated at astrophysical 
sites that are related to baryonic matter, or produced via decay of dark matter 
particles. In both cases the distribution of galaxies (luminous matter) should 
be an acceptable, if certainly not optimal, tracer of the sources. 

Once the previously described scenario is built, test particles can be in- 
jected at the sources and propagated through the intergalactic medium and 
intervening IGMF to the detector at Earth. 

Figure !2Tk-b show the arrival probability distribution of UHECR protons 
as a function of galactic coordinates for a distribution of sources following the 
distribution of luminous matter inside 100 Mpc(CfA2 catalog). A power law 
injection energy spectrum at the sources is assumed, dNi n j/dE oc E , with 
(a) E inj > 4 x 10 19 eV and (b) E inj > 10 20 eV respectively. 

It can be seen that, in contrast to the laminar IGMF case, in this scenario 
information regarding the large scale distribution of the sources inside the 
GZK-sphere can be easily recoverable. The supergalactic plane and the Virgo 
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Fig. 21. Arrival probability distribution of protons (linear scale) as a function 
of galactic coordinates for a distribution of sources following the distribution of 
luminous matter inside 100 Mpc. dNi„j/dE oc E~ 2 , with (a) E in j > 4 x 10 19 eV 
and (b) E inj > 10 20 eV. 

cluster, in particular, are clearly visible between I ~ — 100. It can also be 
appreciated the increase in resolution as the energy reaches the 100 EeV range 
and the gyroradii of UHECR protons become comparable to the size of the 
GZK-sphere. It is also in the cellular model that the deflection angle of the 
incoming particle with respect to the true angular position of the source (see 
figure EH is small enough for an UHECR astronomy to develop at the largest 
energies. 
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Fig. 22. Median and 63% and 95% C.L. for the deflection angle of an incoming 
UHECR proton with respect to the true angular position of the source for the 
example in figure I2l1 All sky average. 



5.3 Anisotropy observations and magnetic fields 

CR anisotropics can be divided, in principle, according to the angular scales 
they affect at a given energy region: (a) large scales, extending across several 
tens of degrees in the sky, (b) medium scales, < 10° and (c) small scales, of 
the order of the angular resolution of the experiment, i.e, ^1—3°. 

At all energies CR are remarkably isotropic. At energies below ~ 100 TeV, 
a first harmonic analysis shows an anisotropy amounting to less than 0.07 %. 
At large energies, on the other hand, measurements are increasingly difficult 
due to lowering statistics and ambiguities in the interpretation of the data due 
to the non- uniformity of the detectors acceptance |106II107| . Nevertheless, all 
the data available is consistent with large scale isotropy (1081 11091 ?.I81|. 

In fact, besides the AG AS A experiment, neither HiRes or Auger have been 
able so far of detecting anisotropy at any energy or angular scale [?. 1811 ITTTj . 

At energies around 1 EeV, i.e., the beginning of the ankle where there 
should still be a sizable Galactic contribution to the flux, Flys Eye has en- 
countered a but statistically significant correlation with the Galactic plane in 
the energy range between 0.2 and 3.2 EeV |112| . They assessed the proba- 



Ultra-high energy cosmic rays: from GeV to ZeV 



31 



bility of this result being a statistical fluctuation of an isotropic distribution 
at < 0.06%. The most significant enhancement is in the interval 0.4-1.0 EeV. 
AGASA, on the other hand has postulated a > 4er excess towards the direction 
of the Galactic center |113| . The excess was confirmed in two independent data 
sets with 18274 and 10933 events in the 1-2 EeV region, with chance prob- 
abilities of 0.3 and 0.5% respectively. The 4.5 a effect observed corresponds 
to 506 events detected in a region of the sky where only 413.6 were expected. 
Associated with this enhancement was a probably dipolar signal towards the 
inner regions of the Galaxy of amplitude 0.04 |114| . An independent confir- 
mation of an anisotropy likely related with this one comes from the SUGAR 
experiment |115| which, different from AGASA, was able to look directly at 
the Galactic center. Unfortunately, these findings have not been confirmed so 
far by either HiRes or Auger |116| . 

At small scales, < 2.5°, comparable with the resolution of the experiment, 
AGASA has postulated the existence of pairs and triplets of events, which 
have grown in number over the years to a total of 7 pairs and 1 triplet (or 
9 pairs if the triplets is counted as 2 pairs) above 40 EeV [?. IHT1 IT23| . 
Since, following AGASA estimations, a total of 1.7 pairs was expected at this 
separation, the results has a chance probability of less than 0.1%. This re- 
sult has not been confirmed by HiRes in the combined AGASA-HiRes data 
set [HI, but still remains a topic of hot debate due to its enormous astro- 
physical significance. It must also be noted that it is difficult to understand 
simultaneously the existence of at least the three original pairs when the actual 
distribution of matter inside the GZK sphere is taken into account |93| . 

Finally, another very promising anisotropic signal coming from the AGASA 
experiment is found as an alignment in the relative orientation of pairs of in- 
coming events above 10 EeV in the Al-Ab plane (Galactic coordinates) on 
scales of < 10° jl 171 1118*] . It must be noted that this anisotropy is funda- 
mentally different from a simple clustering of events in a given angular scale, 
since it is limited to an aligned structure in the two point correlation function. 
This signal can only be produced as the result of charged CR bending their 
trajectories in the Galactic magnetic field. The astrophysical implications of 
this observation have been analyzed in detail in |119j . CR polarization, if con- 
firmed, could turn into a powerful tool for the determination of the number 
of nearby CR point sources and for imposing constraints on the intensity and 
topology of the Galactic and extragalactic magnetic fields. 

Summing up the results on anisotropy at the highest energies so far, and re- 
membering our discussion about the effects of different spacial structures and 
intensities of the IGMF fsections 15 . II and 15 . 2(1 . the high degree of isotropy ob- 
served so far by most experiments, seems to favor a laminar IGMF structure. 
However, it must be remembered that local coherent magnetized structures, 
as our own Halo could be, may be able to de-focus particles coming from point 
sources into an apparently isotropic flux [1251 1126| further complicating the 
analysis. 
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